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ABSTRACT 

We derive the condition on f{R) gravities that admit Kilhng spinor equations and 
construct exphcit such examples. The Kilhng spinor equations can be used to reduce the 
fourth-order differential equations of motion to the first order for both the domain wall 
and FLRW cosmological solutions. We obtain exact "BPS" domain walls that describe the 
smooth Randall-Sundrum II, AdS wormholes and the RG flow from IR to UV. We also 
obtain exact smooth cosmological solutions that describe the evolution from an inflationary 
starting point with a larger cosmological constant to an ever-expanding universe with a 
smaller cosmological constant. In addition, We find exact smooth solutions of pre-big bang 
models, bouncing or crunching universes. An important feature is that the scalar curvature 
R of all these metrics is varying rather than a constant. Another intriguing feature is that 
there are two different f{R) gravities that give rise to the same "BPS" solution. We also 
study linearized f{R) gravities in (A)dS vacua. 
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1 Introduction 



Modifying Einstein gravity with higher-order curvature invariants goes back to early days 
of General Relativity, notably by Eddington and Weyl simply as an exercise of intellectual 
curiosity [1]. A natural and perhaps the simplest generalization of Einstein gravity is to 
replace the Ricci scalar R in the Einstein-Hilbert action with an arbitrary function of R 
[21 [3l m [5]. The resulting f{R) theories of gravity can admit Einstein metrics with a 
positive, negative or zero cosmological constant. This makes f{R) gravities versatile in 
both reproducing Newtonian gravity and studying cosmology, although a convincing f{R) 
theory that fits all the observational data remains elusive. It can be shown [6l[7] that f{R) 
gravity is equivalent to some special class of the Brans-Dicke gravity /scalar theory jS] by 
the Legendre transformation; however, the conversion requires to find the inverse function 
of f'{R) which may not have a close analytical form. Thus in general f{R) gravity should 
be studied on its own right. Since an inflationary model involving the quadratic Ricci scalar 
was constructed first time in [9], there has been ongoing interest in f{R) gravities for the 
last three decades, and the application has been focused largely on the area of cosmology. 
(See, e.g., reviews [101 ttH IE].) 

The fact that AdS spacetimes arise naturally in f{R) gravities suggests that they can 
also be used for investigating the AdS/CFT correspondence [131 (HI [15] . However, apart 
from those with constant scalar curvature, exact solutions in f{R) gravities are difficult to 
come by. When the Maxwell field is introduced in the theory, charged black holes in four 
dimensions were obtained recently [16 ^ \17 \ \TE\. This is possible only because the Maxwell 
field in four dimensions does not contribute to the trace of the Energy-momentum tensor 
so that the Ricci scalar R in the four-dimensional charged black holes remains constant. 
It is easy to see that if f{R) gravity coupled to a matter system whose energy- momentum 
tensor has vanishing trace, the theory admits analogous solutions as those in Einstein gravity 
coupled to the matter. Such solutions with constant R can be viewed as somewhat trivial 
since they do not explore the nature and properties of the function /. Non-trivial solutions 
with varying R have been hitherto unknown in f{R) gravities. With few interesting and 
non-trivial exact solutions, the effort in applying the AdS/CFT correspondence in f{R) 
gravities has been severely limited. 

In supergravities, owing to the existence of Killing spinor equations, much wider classes 
of exact BPS solutions have been constructed. This is because the Killing spinor equations 
can be loosely viewed as the first integral of Einstein's second-order equations of motion, and 
hence they significantly simplify the equations. Killing spinor equations are not exclusive 
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for supergravities. Einstein theory of gravity with or without a cosmological constant in 
any dimension admits a Killing spinor equation^ Killing spinor equations were also known 
to exist in the gravity /scalar system where the scalar potential is constructed in terms of a 
superpotential [19]. It has been recently demonstrated that even when involving form fields, 
some non-supersymmetric gravity theories can admit Killing spinor equations [20^ [2H [22] . 
Examples are rare, but the low-energy effective action of the bosonic string up to the a' 
order was shown to admit Killing spinor equations in any dimension. These Killing spinor 
equations allow one to find new classes of solutions in these theories. 

Obviously, a generic f{R) theory does not admit Killing spinor equations. On the other 
hand, there must exist subclasses of f{R) gravities that do. The simplest example is the 
aforementioned Einstein gravity with/out a cosmological constant. In this paper, we follow 
the technique developed in [201 [21] [22] and propose the two Killing spinor equations for 
f{R) gravities. (Our paper deals only with the f{R) theories in the metric formalism.) 
They involve two functions W and U of R. We then derive the condition on f{R) so that 
the theory admits these equations. We find that for a given choice of W, the function 
/ satisfies a second-order linear differential equation, implying that there exist one (non- 
trivial) parameter family of f{R) gravities for the same W. The function U is then fully 
determined by W and /. Although there are no analytical solutions to the second-order 
differential equation in general, we find many explicit examples of f{R) gravities that do 
admit Killing spinor equations. 

The advantage of having Killing spinor equations in our f{R) gravities is that they can 
be used to construct a large class of solutions involving only first-order differential equations. 
These solutions are analogous to the BPS solutions in supergravities, and hence we refer 
them as "BPS" even though our f{R) gravities are not supersymmetric. The focus of 
our construction is the static domain wall and the FLRW (Friedmann-Lemaitre- Robertson- 
Walker) cosmological solutions, both of which are conformally fiat and of cohomogeneity 
one. It turns out that these "BPS" solutions are solely determined by the function W in the 
Killing spinor equations. It follows that there are two different f{R) gravities that give rise 
to the same solution. We use explicit examples to compare the pros and cons of such two 

^The existence of Killing spinor equations in Einstein gravity or even in supergravity does not imply 
that all background solutions have Killing spinors which are local solutions of the Killing spinor equations. 
Furthermore, manifolds without spin structure would not admit spinors at all. The backgrounds that admit 
Killing spinors are called BPS solutions in supergravities, whilst manifolds with Killing spinors in Einstein 
theory are referred to as ones with the reduced holonomy. The criteria on a theory that admits Killing 
spinor equations will be discussed in section 2.3. 
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f{R) theories. Another characteristic of our f{R) gravities is that they admit at least two 
(A)dS vacua, and the solutions we construct are typically smooth, running from one (A)dS 
vacuum to the other. They provide excellent examples for studying either the AdS/CFT 
correspondence or cosmology. It should be pointed out that although multiple (A)dS vacua 
exist also in Lovelock gravities with the Gauss-Bonnet type of topological terms, there is no 
known flow that links these vacua. This is one major difference between our f{R) gravities 
and the Lovelock type of theories. 

This paper is organized as follows. In section 2, we first give a quick review of f{R) 
gravities and the equations of motion. We then demonstrate a phenomenon with an explicit 
example that there can be new types of (A)dS vacua with vanishing f{Ro) but divergent 
f'{Ro). The same vacuum can be embedded in a different theory with no such singular 
behavior. We then propose Killing spinor equations and derive the T-matrix projected 
integrability conditions. These allow us to derive the condition on / that admits Killing 
spinor equations. We then give a few examples of such theories. More examples will be 
given in subsequent sections. In section 3, we consider "BPS" domain wall solutions. We 
obtain a class of exact solutions that connect two AdS vacua. These solutions include the 
smooth Randall Sundrum II, AdS wormholes and the RG flow from the IR (infrared) region 
to the UV (ultraviolet) region. These explicit solutions demonstrate that f{R) gravities are 
quite attuned to the investigation of the AdS/CFT correspondence. 

In section 4, we examine the FLRW solution with flat spatial directions. We find that 
Killing spinor equations can also be useful to simplify the time-dependent equations. We 
obtain large classes of exact cosmological solutions. One type describes the evolution from 
an inflationary starting point to end with an ever-lasting expanding universe, very much like 
our Universe. We also find smooth metrics of pre-big bang models, bouncing and crunching 
universes. 

In section 5, we give a couple of examples of converting our f{R) gravities to the Brans- 
Dicke theory. In these examples, the inverse function of f'{R) can be obtained as simple 
analytical functions. The majority of our f{R) gravities we obtained in this paper do 
not give rise to a close form scalar potential in the corresponding Brans-Dicke theory. In 
section 6, we study the linear spectrum of f{R) gravities in (A)dS vacua. As one would 
expect, in general the spectrum consists of the massless spin-2 graviton and a massive trace 
scalar mode. We derive the ghost-free and tachyon-free conditions. We discuss the special 
circumstance where the spectrum becomes less straightforward with the kinetic terms for 
the graviton and/or the scalar modes dropped from the linearized action. One intriguing 
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feature of our construction is that there typically exist two f{R) gravities for a given "BPS" 
solution. For such a solution that connects to two different (A)dS vacua, each of the two 
theories is suitable for one of the two vacua respectively. We conclude the paper in section 
7. 

2 f{R) gravities with Killing spinor equations 

2.1 Lagrangian and equations of motion 

The Lagrangian for f{R) gravity in D dimensions is 

Cd = V^f{R) , (1) 

where / is a generic real function. In this paper, we shall be concerned with only f{R) 
theories in the metric formalism, and hence the equations of motion from the variation of 
gf,u are given by 

= F{R)R^, - \f{R)gf,u + - V^V,)F(i?) = , (2) 

where F{R) = f'{R). Note that in this paper, we always use a prime to denote a derivative 
with respect to i?, unless an explicit new variable is given. Taking the trace, we have 

n = RF -\Df + {D -l)nF = {). (3) 

The equations of motion ([2]) can be equivalently expressed as 

n^u = R^u - ^y^^^uF + ^^^^^(/ - 2RF)g^, = . (4) 

Note that although f{R) gravity can be related to the Brans-Dicke theory by the Legendre 
transformation, the resulting scalar potential involves the inverse function of F, which may 
not have a close analytical expression. Thus for general f{R) gravities, the theories are best 
studied in their original forms rather than converting them to the corresponding unnatural 
gravity/scalar system. We shall come back to this point in section 5. 

2.2 New (A)dS vacua in f{R) gravities 

The simplest class of solutions for f{R) gravities are perhaps the metrics with constant R, 
which we denote as Rq. In general, the metrics are Einstein, i.e. R^^ = h.g^y, with the 
effective cosmological constant A = Rq/D. It follows from ([3|) that 

2RoF{Ro) = Df{Ro) . (5) 
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Depending on whether A is positive, 0, or negative, the vacuum solution is de Sitter (dS), 
Minkowski and anti-de Sitter (AdS) respectively. In the special case when both f{Ro) and 
F{Rq) vanish, the equations of motion ([2]) are satisfied simply hy R = Rq. This degenerate 
case allows any metric with constant scalar curvature Rq to be a solution, including some 
Lifshitz black holes |23l I24j . As we shall see in section 6, such an (A)dS vacuum has no 
propagating spin-2 graviton mode, but only a scalar trace mode. (See |25] for a review on 
solutions in f{R) theories.) 

In this subsection, we demonstrate that new classes of (A)dS solutions can emerge in 
f{R) gravity, which are not solutions of but characterized by the divergent F{Rq). To 
illustrate this, let us consider the following f{R) theory 

£4 = (Ti7^i?V48/32 - R , (6) 

It follows from ([5]) that the theory has two vacua with 

Ro = and Rq = 96/3^ . (7) 

It is clear that the two vacua cannot be connected, since for a given cji, the Lagrangian 
can only be real either at the vicinity of i?o = or Rq = 96/32, j-^^^ both. As we 

shall see later, this theory admits Killing spinor equations which enable us to find an exact 
cosmological solution with varying R: 

dsl = —dt'^ + a'^{dx1 + dx2 + dx1), 

It is a straightforward exercise to verify that the metric satisfies ([2]). Let us consider the 
case with /? > 0. The metric is Minkowski when t — )• — 00, and it runs to the de Sitter 
spacetime with Rq = 48/32 when t — )• +00. The Ricci scalar increases monotonically with 
respect to t. The derivation of this solution can be found in section 4.4. 

The solution (j8]) describes a pre-big bang model without singularity. Long before the 
inflation starts, the universe is Minkowski under the perturbative Lagrangian of Einstein 
gravity with higher-order Ricci curvature terms 

C, = 4V3PaV^{R-^ + ---). (9) 

The universe bursts into inflation around t = tQ when the non-perturbative effect takes 
place. 
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What we would like to draw attention to here is that if we rescale Xi — t- e^'^*" Xj and then 
send the integration constant to ~^ — oo, we arrive at a dS metric 

ds^ = -dt^ + e^'^\dxl + dxl + dxl), (10) 

with Rq = 48/3^. It is clear that this is a limiting solution of our f{R) gravity; however, it 
is not a solution of ([5]), which gives only ([7]). This solution is characterized by f{Ro) = 
and F{Rq) = oo, and hence the equation ([3]) breaks down when R ^ Rq. 

In this paper, we are able to uncover such (A)dS solutions because our Killing spinor 
equations allow us to find exact solutions with non-constant R. New (A)dS metrics emerge 
when we let R run to such Rq. It is not obvious to us how to find such a solution in the 
situation when we cannot obtain an exact solution with running R. It is a subject worth 
further investigation. It should be pointed out that from the point of view of the Brans- 
Dicke theory, the solution should be considered as singular since the F{R) corresponds to 
the scalar mode and it blows up at R = Rq. One can also take a different point of view 
that such a theory is intrinsically pure gravity. 

Interestingly, we find that there is another quite different f{R) theory that can give rises 
to exactly the same cosmological solution ([8]), namely 



£4 = a2^R 12/3 - x/3(48/32 - R) arctanh^^ "^^^^^ ^ ) . (11) 

This somewhat more complicated theory is convergent at Rq = 48/3^, namely 

F(i?o) = 24/3^2 , Rof'iRo) = 16/3(72 . (12) 

Thus in this case, the equation ^ is satisfied. As we shall discuss in section 6, this theory 
is ghost free in this de Sitter vacuum. On the other hand, the theory becomes singular at 
i? = 0, since f{R) ~ RlogR when — )• 0. Thus, F is divergent at i? = in this case. 

Thus we find an interesting phenomenon in f{R) gravities. The same cosmology ([8]) can 
be generated by two different classical actions, but neither theory can smoothly describe 
the full evolution. In the perturbative flat region, the Lagrangian ([6]) is a better theory. 
The later inflationary epoch is better studied by the theory (jlip . This is similar to the 
phenomenon in differential geometry that a manifold typically requires multiple different 
but overlapping coordinate patches to cover it. More detailed analysis will be given in 
section 6. 
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2.3 Killing spinor equations 

In the previous subsection, we present two f{R) gravities that give rise to the same cosmo- 
logical solution with varying R. Such an exact solution is possible because the two theories 
are special in that they admit Killing spinor equations. In this subsection, we derive the 
condition on / so that the f{R) theories admit Killing spinor equations. 

It is clear that for a generic function / there can be no consistent Killing spinor equations. 
However, for some classes of functions, the theories can admit Killing spinor equations. The 
simplest example is that / = R— {D — 2)Aq, whose Killing spinor equation is well established 
and given by 

D,e^[D, + l^/^^T,)e = 0, (13) 
where is a covariant derivative on a spinor, defined by 

D^e = d^e + lco';^'Tabe. (14) 

Whilst the equation D^e = has a clear geometric interpretation that e is a covariant 
constant spinor, the extra term in (I13p lacks an immediate explanation. One defin- 
ing property for the Killing spinor equation is that the F-matrix projected integrability 
condition gives rise to the Einstein equations of motion, namely 

= r^[^^, D,]e = ^T^'iR^, - Aog^,)e . (15) 

Thus we see that the projected integrability condition is satisfied provided that the metric is 
Einstein with cosmological constant Aq. It follows that without the extra term in (|13p . 
the Killing spinor equation would become irrelevant to the theory with the cosmological 
constant. Conversely, the significance of the above projected integrability condition is the 
following. Although the existence of Killing spinor (jl3p for a background does not in general 
imply that it always satisfies the equations of motion R^^, = Aqq^^, , if one of the Killing vec- 
tors constructed from the Killing spinors is time-like, the background then indeed satisfies. 
This provides a powerful tool for constructing exact solutions, since first-order equations are 
much more manageable than the second-order ones. (The integrability condition without 
the F-matrix projection is related to Riemann tensor, and hence not related directly to the 
Einstein equations of motion. This implies that not all solutions have Killing spinors, even 
if the theory admits the Killing spinor equation.) 

We are now in the position to derive the condition on / so that the theory admits Killing 
spinor equations, whose F-matrix projected integrability conditions are analogous to (|15p . 
Such Killing spinor equations were constructed recently for non-supersymmetric theories 
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involving the metric, a dilaton and a form field [20\ [2H [22] . It turns out the condition is 
very restrictive and the only known non-trivial examples are either the low-energy effective 
action of the bosonic string up to the a' order |20[ [2T] or the Kaluza-Klein theory with a 
non-trivial scalar potential [22]. In f{R) gravity, we propose 

V^e= (^D^ + W{F{R))r^y = 0, (r^V^F + C/(F(i?))) e = , (16) 

where the functions W and U are to be determined. The first equation is the natural 
generalization of ()13p . The second equation is inspired by the fact that f{R) gravities are 
effectively a special, albeit inconvenient, class of Brans-Dicke theories. In supergravities, 
the left-hand sides of equations in (jl6p would be the supersymmetric variations for the 
gravitino and dilatino fields, and F, W are analogous to the scalar and the superpotential 
for the scalar. If / is linear in R and hence is a constant, we must set U = and we 
recover the previous example. 

To establish the relevance of these two Killing spinor equations (jl6p to f{R) gravities, 
we follow the procedure developed in [201 [211 [22] . We first act on the second equation with 
r'^Vj/, which gives 

[DF -U{U + 2{D -l)W))e = 0. (17) 

In this paper, a dot is always denoted as a derivative with respect to F. Note that for any 
function X, we have X = X' /f" , by the virtue of the Leibnitz' chain rule. The F-matrix 
projected integrability condition for the first equation in ()16p is somewhat more involved, 
and we find it is given by 

Tf'iv^Me = i^'^7^^,.e + ^r^x^. , (18) 

where TZ^u is given by (jl]) and X^jy is given by 



X, 



,2 4(L» - 2)WWU 1 



2WU + 4(1. - l)W^ + - 2^^) j^- 

yF 2W -U ) ^ ^ ' 

Thus for an f{R) theory to admit consistent Killing spinor equations (|16p . we must have 

21^ - if ■ 1 

U{U + 2{D-l)W) = —^—-{Df-2RF), 



2{D-2)W' ' ' ' ' 2(L»-1) 

2WU + A{D - l)W^ + ^ - ^^^1^(/ - 2RF) = . (20) 

These equations can be reduced, giving rise to 

AD{D-l)W^ + R ■ _ {AD{D-l)W^ + R)W 

^{D-l)W ' {A{D-l){D-2)W^ + R)F-f ^ ' 
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Thus, once W and / are known, U comes out straightforwardly. It is advantageous to 
express all functions in terms of the variable R, in which case, the second equation above 
becomes 



One way to view this equation is that it is a non-linear differential equation of W for 
a given /. Owing to the non-linearity, however, a solution for W does not always exist 
for a generic function /. An alternatively view point is that (j22p is a second-order linear 
differential equation for /. Thus (I22p must have solutions of / for any given W, even 
though the analytical expression for / may not exist. This is very much parallel to the 
gravity /scalar system, where a generic scalar potential may not admit a superpotential, 
but any superpotential can yield a potential. The difference is that here there is no simple 
expression directly for / in terms of function W, but instead / has to be solved via the 
second-order linear differential equation. What is curious is that for a given W, there can 
exist a two-parameter (one non-trivial) family of f{R) theories, since second-order linear 
differential equations tend to give two solutions associated with two integration constants. 
In other words, there are two different f{R) gravities for one W. Of course the function U is 
also different. This is different from the usual gravity /scalar system whose scalar potential 
can be expressed simply in terms of a superpotential when the theory admits Killing spinor 
equations. In section 2.2, we have demonstrated such two different f{R) gravities can give 
rise to the same solution. Furthermore, the cosmology at two different evolution epochs is 
better studied by either one or the other f{R) theories. 

To summarize, we give a simple way of constructing f{R) gravity with Killing spinor 
equations. We can begin with a function W, from which the f{R) is determined by (j22p . 
and the expression for U follows straightforwardly from (j2ip . The corresponding Killing 
spinor equations are then (jl6p . 

It should be emphasized that the Killing spinor equations we have discussed are for the 
theory with generic backgrounds rather than for a specific solution. For an f{R) theory 
with such Killing spinor equations, we can check whether a specific solution preserves Killing 
spinors and obtain the fraction of the maximally allowed Killing spinors that survive. For 
example, we can establish whether the (A)dS vacua satisfying ([5]) admit any Killing spinor 
that satisfies ()16p . Furthermore, armed with these Killing spinor equations, we can construct 
new solutions. As mentioned earlier, the Killing spinor equations (|16p can be viewed as the 
first integrals of the Einstein equations of motion. For a background with a Killing spinor 




(22) 
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e satisfying (jl6p . the integrability conditions imply, as we have shown, 

ne = 0, T^'n^.^e = , (23) 

where TZ and TZ^u are defined in ([3]) and (jl|) respectively. Thus we see immediately that 
the trace equation ([3]) is automatically satisfied, and the second equation in ()23p takes the 
same form as (|15p . Note that the first equation above is implied by the second equation, 
and hence it is not independent. If the Killing vector constructed from the Killing spinor 
e, namely K^{e) = eF^e, is e.g. time-like, it is straightforward to demonstrate that the 
equations of motion are satisfied as well. Thus for our f{R) gravity with Killing spinor 
equations, we can obtain a class of solutions by solving only the Killing spinor equations 
(jl6p with the corresponding time-like Killing vector, instead of solving the more difficult 
Einstein equations directly. This class of solutions are analogous to the BPS solutions in 
supergravities, and we shall refer them as "BPS" even though our f{R) gravities are not 
supersymmetric. Note that in practice, our task may not be about constructing the most 
general "BPS" solutions, but only some special solutions with some simple ansatz. In this 
case, the existence of some Killing spinors that satisfy the Killing spinor equations can 
determine the ansatz fully. We can then simply substitute the result into ([2]) to verify 
whether it is a solution or not. 

2.4 A few examples 

Here we give some simple explicit examples of f{R) gravities that admit Killing spinor 
equations (jl6p and give the corresponding W and U functions. An obvious search is to 
consider quadratic f{R) gravity in D dimensions, namely 

f{R) = (tR-{D- 2)Ao + aR^ , (24) 

where o", a and Aq are constants. For these general parameters, we find that the Killing 
spinor equations do not exist. However, if the three parameters satisfy the following con- 
straint 

4(D - 1)2^2 

"^° = -(2)-2)(5I)-2)2' ^2^) 
we find that they do exist, with 

U = -^^(D{D + 2)a + \{D-A){bD-2)F)F\, 
^ = -l^{iD + 2)a-{f>D-2)F)F-'2^ 



2 



3 



IQ{D + 2){D - l)a ■ 
12 



(26) 



Another example is also a £)-dimensional theory: 

f{R) = aR{R + p)^ . (27) 
We find that the corresponding W and U are given by 



o {D-2)R{{D-1)R + I3D] 
W = , , U = ^ . (28) 

In section 2.2, we have shown that the two four-dimensional f{R) gravities ([6]) and 
(jlip give rise to the same cosmological solution Both theories admit Killing spinor 
equations, with the same i.e. 

W = ^. (29) 



However, these two theories have diff'erent U functions. For the theory ([6]), the function U 
is 

16/3^48/32 - R 



For the theory (jlip . we have 

£12(12/3(32/32 - R) ^48;S2 _ ^ + ^ i?(64/32 - i;)arccothf 

U = ^ , Wl!^-R^_ 

16/3v^48/32 - ii 

More examples of f{R) gravities that admit Killing spinor equations will be given in sub- 
sequent sections. 



3 "BPS" domain wall solutions 
3.1 General properties 

As discussed in section 2, Einstein metrics arise naturally in f{R) gravities. However, for 
a generic f{R) theory, any exact solution beyond constant R is more or less impossible to 
construct, since one has to handle in general fourth-order non- linear differential equations. 
However, for our f{R) gravities with Killing spinor equations, exact solutions with varying 
R can be obtained. In particular, we shall consider static domain wall solutions, with the 
ansatz 

dsl = dr^ + e'^^^''^dx^'dx^. (32) 

Making a natural choice of vielbein, e"^ = dr, e* = e^dx'', we find that the non- vanishing 
spin connection is = A^^e^ . The Ricci tensor and scalar are given by 

= -(£'-l)(Arr + ^',), R,., = -{A^„ ^ {D - \)A\)g^, , 
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R = -2{D - I) A ^rr - D{D - I) A^^ . (33) 



For the metric ansatz ()32p . the KiUing spinor equations (|16p can be easily solved. The 
Killing spinors are given by 

e = es^eo , T^eo = eo , (34) 

where eq is a constant spinor in the {D — l)-dimensional world-volume. The full set of 
equations of motion is now reduced to simply 

A^r = -2W . (35) 

This is a tremendous simplification of the Einstein equation We see that the struc- 
ture of domain wall is completely determined by the function W. This is analogous to 
supergravities, where the BPS domain walls are solely determined by the superpotential. 

From the Killing spinors (j34p . we can construct the Killing vectors, namely K^^ (eq) = 
e^eoF^^eo, and hence K"" = e^eoeo, = e^eoF*eo and i^* = e^eoF'eo. Whether this Killing 
vector can be time-like or not depends on dimensions and the F-matrix properties. Since eg 
can be any constant spinor in the D — 1 dimensional spacetime dx^dx^^ we would expect 
that we could choose some appropriate eo such that = = K^, but 7^ 0. In the case 
when D is odd, and hence [D — 1) is even, such an eo can be easily found. For example, 
let us consider the convention that e = e^F* with Hermitian F* and anti-Hermitian F*, the 
Killing spinors eg satisfying the above conditions are given by 

eo = (1 + i)m , (36) 

where 7 = Y\ - F* with a so chosen that 7^ = 1. It is straightforward to verify that 
the corresponding if is a time-like Killing vector. Following the discussion in the previous 
section, the equation (|35p must satisfy In even D dimensions, the discussion is somewhat 
more complicated. Let us consider two constant spinors Cq = (1 it Fi)?70) and denote 
as the corresponding Killing vectors. It is clear that are both null vectors and hence 
we can choose a convention such that K = + K~ is time-like. Since we must have 
TZn^K'^ = 0, it is then straightforward to show that TZf^^, = 0. Note that this demonstration 
works in odd D dimensions as well. 

In fact, for such a simple background, it is quite easy to demonstrate that (|35p indeed 
satisfies all the equations of motion by simply substituting ()35p directly into ([2|). From the 
expression for the Ricci scalar in (133p . we obtain the first-order equation 
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Together with (p2|) . we find that 

p _ {R + ^D-l){D-2)W^)F-f 

A{D-l)V ■ ^ ' 

It is now straightforward to estabhsh that the full set of equations of motion are all satisfied. 

Now let us consider the general properties of the solution. We use X to denote the 
right-hand side of (f37l) . namely 



Here we treat i? as a function of W. The equation (|37p is of the first order, and can be 
solved as 

If X{W) has a zero, such that 

X{W) = X'{Wo)iW -Wq) + --- , (41) 
we find that near the region of Wq, the solution is given by 

, X'(Wo)r 

W = Wo + e''''^'^^^\ e^ = e^p[-Wor-^^^j. (42) 

Thus if X'(Wq) > 0, the metric becomes AdS when r — ?• — oo. The resulting metric is AdS 
horizon if Wq < 0, and it is AdS asymptotic boundary if Wq > 0. If on the other hand we 
have X'{Wq) < 0, the metric becomes AdS when r — )• +oo. The resulting metric is AdS 
horizon if Wq > 0, and it is AdS asymptotic boundary if Wq < 0. Thus we see that near 
the region oi W = Wq, the solution is regular, approaching either the AdS horizon or the 
AdS boundary. Thus, If X{W) has at least two roots, we can expect smooth solutions that 
run from one AdS to the other, associated with two adjacent roots. 

3.2 A class of exact solutions 

As discussed above, in order to construct smooth solutions, it is necessary that the function 
X(W) has two roots. In this subsection, we consider a class of W, which is given by 



R = 4:D{D-l)i^{a-l)W^ + bW + cj. (43) 
where a,b,c are constants. It follows that X(W) is quadratic: 

X{W) = D{aW^ + bW + c), (44) 
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which has two roots, given by 



2W^ = A± = . (45) 

a 



We require that the discriminant A = 6^ — 4ac > so that the two roots are real. The W 
equation (|37l) imphes that 

= -— (&+ \/A tanh(iL»^(r + ro))) . (46) 



It follows from (j35p that we have an explicit solution 

g2A ^ (^^bDr cosh2(lL)^/A (r - ro))) . (47) 

The metric is smooth with r running from — oo to +oo. In both limits, the metric approaches 
AdS, namely 

g2A^g2A±r^ for r^±oo. (48) 

The resulting effective cosmological constants in the AdS limits are A.± = —{D — 1)A^. 

For the choice of W given by (j43p , it follows from ()22p that / can be determined by the 
following differential equation 

D{aW^ + hW + c)(2(a - l)W + b)Wf,ww 

-(A{a - l)W{{aD - 1)W'^ + cD) + b{{5aD -AD- 2)W'^ + cD) + b^DW^f^w 
+D{2{a-l)W + hff = Q. (49) 

For generic parameters (a,6, c), there is no analytical solution for /; however, for some 
special choices of these constants, we obtain explicit f{R) gravities. One example is that 
6 = 0. In this case, the f{R) is given by two hyper geometric functions 

f = aiW' 2Fi(x-, x+; --; -^) + a2 2Fi{y.,y+; ^; -^) , (50) 

2 c 2 c 

where o"i and o"2 are integration constants properly chosen so that the function is real, and 



2 - 3aD ± JA + aD((a + 8)D - 12) 

= 4aD 

2 + 3aD±JA + aD{{a + 8)D - 12) 

= 4^ • 



The expression for W in terms of R can be obtained from (143 h . Another example is provided 
that a = 1 and c = 0, in which case, we have 

f{R) = aR(R + AD{D -l)b^) . (52) 



Note that here we have not presented the other choice for / which involves hyper geometric 
functions. This simple f{R) theory were presented earlier in section [27il 
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3.3 Randall-Sundrum II 

The Randall-Sundrum (RS) II scenario is characterized by the metric profile of the type 
e^"^ = e~^'^l''l, with positive constant k [26]. In other words, the metric approaches the 
AdS horizons at both r — )■ ±00 limits, with one maximum in the middle. (Although r is 
a non-compact coordinate, the volume integration is finite with respect to r.) This can be 
achieved in our domain wall solution by requiring a < 0, c < 0, which ensures that 

A+ < , A_ > . (53) 

Note that the constant b can be arbitrary in this case. Our solution describes in general the 
asymmetric RS II scenario when A+ -|- A_ does not vanish, since the cosmological constants 
A± are not equal at the two AdS horizons. This of course is not crucial for trapping gravity 
on the wall. A symmetric RS II can be obtained by further requiring 6 = 0, whose f{R) 
theory is given by ([50]) . 

To study the trapping of gravity on the wall, it is advantageous to express first the 
metric in the conformally-flat frame, namely 

= e2^(^) {dx^^dx^ + dz^), (54) 

where the coordinate z is related to r as follows 



V cosh(iD^/Ar) / aD^/A aD' aD^/A 

Note that we have set the inessential tq to zero. The linear fluctuation of the graviton modes 
in the {D — l)-dimensional flat world- volume in the context of f{R) theory has not been 
studied yet except for D = 5 [27]. Thus we shall focus our attention on five dimensions. 
Following the procedure outlined in [27], we let h^i, = e~^'^/^F~^/^n^,^(/)(z), where n^j^ is 
transverse and traceless, we find that 

-l^,zz + Vct> = 0, (56) 

where the Schrodinger potential F = Vq -|- Vi contains two parts. The first part is given by 

^0 = + liD - 2)A,,, + 1{D - 2fA% , (57) 



which is the same as that in [26]. For our general domain wall solution, we find 

D-2 

exp ^^[t 

, 2 



^0 = ^k^ + ^^exp 



l{hr+^ log cosh(iL>VA (r - ro))^ 
(^aZ)(^sech(iD\/Ar))^ + (^6 + VAtanh(iL>VAr))) . (58) 
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The second part is the contribution from F{R), and it was estabhshed only for D = 5 |27j . 
It is given by 

Note that the conversion to the coordinate z does not give a close form for general pa- 
rameters and one can appeal to the numerical approach. For some choice of parameters, 
the explicit y as a function of z can be obtained. For example, let a = —2/D,h = and 
c = D/2, in which case, we have 

cosh(L'r) 

It follows that r = 4arcsinh(£'z), and hence 



1,2 , D^ {D-2){D^z^- 2) 
SiD^z"^ + 1)2 



Vo = kk' + 'o^^T, ;r ' ■ (61) 



This potential profile is very much like the one obtained in [29] and is capable of trapping 
gravity on the wall, but ours is realized by f{R) gravity, with 

fiR) = ai|l^|SFi(l(l-^/2^), 1(1 + ^/2^);!;-^) 

, , , , , AW"^ 

+CJ2 2Fi(-l - ^\/2^, -1 + i\/2^; -1; -^) , 

^ 2DHD-1)-R 

4(Z)- l)(D + 2) ' ^ ^ 

The alarming-looking complex arguments in the hypergeometric functions do not prevent 
the f{R) from being real provided that \ W\ < D/2. Of course, in order to demonstrate the 
trapping of gravity in the f{R) theory, we also need to look at the contribution from Vi. 
Let us consider the simpler example associated with a2 in (I62p in D = 5. We find that 

25 



Vi 



(1 + 25^2)2 



6 - 3tanh2 (y/s arcsm 



Vl + 25z2 

_ /3-h75z2 / ^ 5z 

-oOz\ ^ tanh v 3 arcsm- 



(63) 



1 + 25^2 V Vl + 25z2- 

The Schrodinger potential V = Vq + Vi becomes a bit more complicated. It has a local 
maximum V = at z = 0, and two negatives minimums when we increase \z\, and it 
becomes positive until it hit a maximum before it approaches zero at \z\ = oo. Thus we 
demonstrate that f{R) gravity can easily reproduce the smooth RS II scenario with thick 
domain walls. Note that we have added in an absolute value symbol on W"^ so that the 
f{R) is a symmetric function of W. This does not create a discontinuity in ()49p at W = 
since it involves only up to the second-order derivatives. The Ricci scalar for these solutions 
runs from Rq = —D^{D — 1) on the AdS horizon at r = — oo to the maximum value of 
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—2D^{D — 1), and then decreases and approaches Rq again on the AdS horizon at r = +00, 
and correspondingly W runs from —D/2 to D/2. It is easy to verify that for (j62p . the 
constant Rq = —D^{D — 1) solution indeed satisfy the equation ([5]). Furthermore, we have 
, , 3L>sinh(iV^D^7r) coshfi^/F^vr) ^ ^ 

The /"(i?o) is divergent unless the constants ui and (T2 are specifically related. Since the 
expression becomes very complicated, we shall only give the result for D = 5. In this case, 
we have 



172 = — — — coth(i\/3 vr) o"! , (65) 



125\/3 
"56 

and 

F{Ro) = ^csch(l V3 ^) , Rof'iRo) = -^csch(i^/3 vr) en , (66) 
where Rq = —500. 

More generally, the RS II scenario arises in any X(VF) given in (|39p with a profile that 
it has one positive, one negative root and a maximum in between. 

Note that an analytical domain- wall solution in quadratic f{R) theory together with 
additional scalar with cj)'^ potential was obtained in [28]. 

3.4 AdS wormholes 

The situation is quite different if we have 

A+ > , A_ < . (67) 

This can be achieved by requiring a > and c < 0, while b can be arbitrary. In this 
case, the solutions describe smooth AdS wormholes that connect two AdS boundaries at 
r — >• ±00, with no bulk singularity and horizon in between. For non-vanishing b, the 
wormhole connects two asymmetric AdS boundaries with different A±. For 6 = 0, the AdS 
boundaries are symmetric. Let us present a relative simple example with a = 2/D,b = 
and c = —D/2. In this case the theory is given by 

a2x 



/ = ai\Wf2F,{-^{l + VD^),-^{l-VD^);-^;-^W' 
+CJ2 2^1(1 - ^VD^, 1 + ^VD^; f ; ^W^) 



^ _2D^2iD-l) + R 

A{D-l){D-2) ^ ' 

We add an absolute sign in W for the same reason explained in the previous subsection. 

The solution is quite simple: 

ds^ = cosh^ {Dr)dx''dXf, + dr'^ . (69) 
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The metric reaches AdS boundaries with Rq = —D^{D — 1). For D = 1 + A;^, the function 
/ is given by simple functions. For example, when = 5, we have 

f[R) = aiy/-{R + 200)3 + o-2(i? + 50)VR + 500 . (70) 



For the f{R) = ai^/-{R + 200f theory, we have F{Ro) = -15\/3o-i, and hence the ([5]) 
is satisfied. For the f{R) = (T2{R + 50)\/R + 500 theory, on the other hand, the F{Rq) is 
divergent. Such a situation was discussed in section 2.2. 

It is worth commenting that although there is a no-go theorem in the usual Einstein 
theory that the configuration with two AdS boundaries connected in the bulk without a 
horizon separating them will violet the energy condition [33]. This no-go theorem can be 
easily circumvented in higher-order derivative theories. Smooth wormholes with two AdS 
boundaries were constructed in Einstein gravity with the Gauss-Bonnet term [30\ \3T\ [32] . 
Our examples demonstrate that wormholes arise naturally in f{R) gravities as well. More 
generally, wormhole solutions occur in any X(W), given in (j39p . with a profile that it has 
one positive and one negative roots with a minimum in between. Note that our wormhole 
solutions are brane-like and static. Only stationary brane-like wormholes were known to 
exist in Einstein gravity and supergravities in higher dimensions \34: \ \35 \ [36] 137]. 

We can also consider a different parametrization. We set a = 1 and c = 0, but with 
non- vanishing b. In this case, one / is given in (|52p . The domain wall solution is given by 



e2^= (l + e^^(^-^°))^. (71) 



For 6 > 0, this solution describes a wormhole that connects a flat spacetime at r — t- — oo to 
the AdS boundary at r — )• +oo with Rq = —AD{D— 1)6^. Note that the solution with 6 < 
is equivalent to 6 > 0, by reversing the sign of r. Similar solutions that connect the AdS 
boundary in one asymptotic region to the flat spacetime in another have also been found 
in supergravities [35 [ I36 [ [37]: however, these solutions are stationary rather than static. 

As discussed earlier, there is a different f{R) theory that would give rise to the same 
solution (j7ip . It is much more complicated, given by 



where W = R/{4:D{D — 1)6). In the case of D = 4, the expression is simpler, given by 

(73) 

This theory satisfies that F(i?o) = 245a and Rof"{Ro) = 165ct. 



f{R) = aR^Ub - VsCRTIsP) arctanh(^^^^^) 
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3.5 RG flow from IR to UV 

In this case, we have A+A_ > 0, which can be achieved by requiring }p' > Aac > 0. It is 
clear that X± being both positive is equivalent to the case with both being negative, by 
merely reversing the sign of the coordinate r. Let us thus discuss the case with both being 
positive. If a > and hence 6 > 0, we then have 

A+ > A- > . (74) 

The metric describes a flow running from the AdS horizon with A_ at r — )■ — oo to the AdS 
boundary with A_|_ at r — )• +oo, which corresponds to the IR and the UV regions in the 
dual conformal field theory respectively. The cosmological constant A_ in the IR region is 
smaller than the A+ in the UV region. This type of behavior is similar to the domain wall 
solutions in supergravities [38] ■ Note that such solutions can be obtained in any X(VF) that 
has two adjacent positive roots with a minimum in between. If instead, a < and hence 
6 < 0, the cosmological constant in the IR region is bigger than that in the UV region. Such 
a solution occurs in any X(W) that has two adjacent positive roots with a maximum in 
between. 

3.6 On holographic c-theorems 

As we see in our explicit constructions of domain wall solutions, f{R) gravities are quite 
suitable for investigating the AdS/CFT correspondence. One natural question is to examine 
the holographic c-theorem. One may view that f{R) gravity is simply Einstein gravity with 
an effective energy-momentum tensor built from the Ricci scalar, with the equations of 
motion ([2]) expressed as 

Gfiu = T^u ~ Gfiv — Qp,v 1 (75) 

where = R^^u — \gji,yR is the Einstein tensor. If one takes this point of view, one can 
follow the procedure in [38] and define 



vr 



(D-l)/2 



air) = — ^ . (76) 

Its variation with respect to the co-moving coordinate r is given by 

„(D-l)/2 

a , = — ( (T<=ff)^ - [T'^^Yr] . (77) 

The holographic c-theorem follows provided that a^r ^ 0, which implies that the cosmologi- 
cal constant at the IR is smaller than that in the UV. It is clear that whether the c-theorem 
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holds or not depends on the specific choice of f{R) gravities, and we have examples that 
both support and violet the c-theorem. 

A different approach is to treat f{R) gravities as Brans-Dicke theory, and the c-theorem 
is then dictated by the corresponding scalar potential. The third approach is treating f{R) 
theory as a pure gravity theory that can coupled to additional matter so that the equations 
are now given by 

Q,.^ = ■ (78) 
This follows the same approach of [39] where all ghost-free curvature squared and cubic 
terms were considered. A monotonic function a can be found in these higher-order theories 
|39| . It is of great interest to investigate the constraints on / so that the holographic 
c-theorem also holds and whether such constraints are consistent with the conditions for 
Killing spinor equations. 



4 "BPS" cosmology 
4.1 The set up 

It is well-known that the de Sitter spacetimes also admit Killing spinors, even though 
Einstein gravity with a positive cosmological constant cannot be supersymmetrized. The 
Killing spinor equation is given by 



D,e^[D, + ^^^^^T,y = Q, (79) 

This property of de Sitter space was exploited in constructing de Sitter "supergravities" 
|401l41j which are effectively the analytical continuation of AdS supergravities. The function 
W in Killing spinor equations (I16p is pure imaginary in this case. We would like to assume 
implicitly that W and U in Killing spinor equations are real. Thus for the purpose of 
studying cosmology, we would like to rewrite the Killing spinor equations as follows 

V^e={D^ + \WT^y = {), (r^V^F + iC/)e = 0, (80) 

where U is given by 

A{D-l)W 

It is important to note that our procedure of sending W and U to imaginary values does 
not affect the reality of the function /, which now satisfies 



r-^^ +1 ^/ = 0. (82) 

[R - AD{D - 1)WAW (R-4:D{D -1)WAW 
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We now construct "BPS" cosmological solutions that admit Killing spinors. The ansatz 
is the FLRW metric with flat spatial directions 

ds'^ = -dt^ + a^dx'dx' . (83) 

Requiring that the solution admit Killing spinors, the full set of Einstein equations of motion 
is reduced to 

^ = 2W , (84) 

a 

which implies the following first-order equation 

«zi)2^. ,35, 

We verify that (j84p indeed satisfies ([2]). Smooth cosmology emerges when y(M^) has two 
adjacent roots corresponding to two de Sitter spaces. The cosmological evolution runs from 
one de Sitter to the other. 

As an illustrative example, let us consider 

R = 4.D{D-l){^{a + l)W'^ + pW + , (86) 

such that Y{W) = D{aW^ + /3W + 7). The corresponding f{R) can be determined by the 
following second-order linear differential equation 

4(a + l){aD + 1)W^ + /3(5aL» + 4D + 2)W^ + L>(4(a + 1)7 + f3'^)W + /3L>7 



DW{2{a + l)W + /3){aW^ + + 7) 



w 



2(a + l)l^ + /3 

+^'^^+W(aH^2 + ^^ + ^)/-0- (87) 
When /3 = 0, the equation can be solved explicitly, giving 

f = aiW' 2Fi{i.,x+; -; ) + ^2 2^1 (y-, y+; -^) • (88) 

27 7 

where ai and (T2 are integration constants properly chosen so that the function is real, and 

3aD - 2 ± ^/4 + aD{{a-8)D + 12) 



x± 



4aD 



-3aD - 2 ± JA + aD{{a-8)D + 12) 

= ^ ■ (89) 

The general "BPS" cosmological solution for (j86p is given by 

a = (e^^'cosh\lDVA {t - to))) ~^ , (90) 

where A = — 4a7 > 0. The solution approaches de Sitter spaces in both t — )■ ±00 limits, 
with a ~ e'^**, where 

A. = -^^. (91) 
a 
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The corresponding cosmological constants of the de Sitter spaces in these hmits are A± = 
{D — 1)A^. Depending on the sign and values of X±, various cosmological scenarios emerge. 

If we take a view that f{R) gravity is simply Einstein gravity with an effective energy- 
momentum tensor, as in (j75p . the effective energy density and pressure for our cosmology 
are given by 



cff _ 



2 

^2 { l3 + ^/Atanh{^ D\^ty 



a? o? 



p 

^ \{aD{D -2) + 2{D -l))A(sech{^DVAt 



2{D - l)a2 

-2(Z)- 1)(/?2 + A + 2/3\/A tanh(iZ)VAt))j . (92) 



Then we have 



^ 2(Z) - 1)1/3 cosh(iD\/At) + \/Asinh(iD\/At) 



Note that in this discussion, we have set the inessential Iq to zero. In the limits of t — )• ±00, 

we have «;°^(iboo) = — 1 as one would have expected. Since A > 0, the sign choice of 

^^eff _|_ 2^ ii^Q "dark energy" depends solely on the parameter a. An extremum occurs 
at 

t = ^^arctanhf^^^^ , (94) 

D^/A V /3 y ' ^ ^ 

corresponding to 



{W^^ + l)extremum = _ ^^j^ • (95) 

4.2 From inflation to ever-expanding universe 

The parameters a, (3 and 7 should be chosen such that A4. and A_ are both positive. Fur- 
thermore we must have A+ < A_. In this model, the the universe starts an inflation with 
a bigger cosmological constant A_ at t — t- —00 and end with an ever-expanding de Sitter 
universe with a smaller cosmological constant A_|_. Since we have 

A+-A_ = -^. (96) 
a 

It follows that we must have a > 0. This implies that for this model the sign choice of 
{w°^ + 1) for the dark energy is always positive throughout the evolution. Furthermore, 
(3 must be negative and 7 must be positive. As a semi-realistic model of our universe, we 
require 

(97) 



A_ 
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This can be achieved by requiring 4a7//3 << 1, in which case we have 

A_~-?^, X,~-^. (98) 

a p 

It can be shown that in the later part of the evolution, f{R) ~ i? — B? / + 0{B?). Note 
that this type of semi-reahstic solutions emerge as long as has two positive roots 

with a positive maximum in between. 

If instead we have A+ > A_ > 0, which can be achieved by requiring q < 0, /3 > and 
7 < 0, the universe would start with a mild inflation, and inflates faster and faster. Such 
a solution arises in general when y(VF) has two positive roots with a negative minimum in 
between. Such a theory provides a model for the multi-stage inflationary scenario. 

4.3 Bouncing universe 

The universe bounces when A_ < 0, but A+ > 0. This occurs when a < and 7 > 0. The 
minimum a occurs when t = tmim given by 

2 /3 

imin — ^0 = ^arctanhf^^) . (99) 

D^fK Va 



The minimum scale factor is given by 



l-^)^(^)^. (100) 



A/ 

This type of bouncing universe emerges when y(VF) has one positive and one negative 
adjacent roots with a positive maximum in between. 



4.4 Pre-big bang model 

In the special limit, namely 7 = and a = — 1, the solution is simple, given by 

a= (^l + e^^(*-*"))*. (101) 

In this case, one f{R) gravity takes a simple form, i.e. 

f{R) = r[aD{D - l)/32 - i?^ ^ . (102) 

We shall not present the other f{R) gravity that gives rise to the exact same solution. 
The solution connects the flat i? = region to the R = 4:D{D — l)/3^ de Sitter space. It 
can be used to model the singularity-free inflation scenarios. In particular it predicts a 
"pre-big bang" flat universe which bursts into inflation by the non-perturbative effect of 
the higher-order curvatures. The four-dimensional case was discussed in section 2. 
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4.5 Smooth crunching universe 

If we have A_ > 0, but A+ < 0, the universe starts with an inflation, but end with a big 
crunch. What is interesting is that usuahy such model encounters a curvature singularity 
at the crunch. But in our solution, the universe shrinks in the manner of a de Sitter space, 
and hence there is no singularity. 

5 Relating to the Brans-Dicke theory 

It is well-known that f{R) gravity can be cast into the form of Brans-Dicke theory by the 
Legendre transformation. To see this, one starts with the Lagrangian 

= V^{f{x) + f,xix)iR - x)) . (103) 

Variation with respect to x gives rise to 

1^^{R-X) = 0. (104) 

Thus provided that /.^^ / 0, it follows that x = R-, a^id hence (jl03p is the usual f{R) 
theory. Alternatively, one can deflne 

^ = /,x(x), (105) 

and hence the f{R) gravity is equivalent to the Brans-Dicke theory of the type 

C = ^g{ipR + fixm-^xiv)) . (106) 

This is a special class of Brans-Dicke theory with no manifest kinetic term for 93. The 
conversion of f{R) gravity to the Brans-Dicke theory requires finding the inverse function 
of F = /', which in general does not have explicit analytical form. In most of our examples 
that admit Killing spinor equations discussed in this paper, the f{R) theories are better 
discussed on their original form, rather than converting to the corresponding Brans-Dicke 
theories. There are couple of examples we find that can be converted into the gravity/scalar 
system, where the scalar potentials are expressed in terms of simple functions. 

5.1 A quadratic f{R) theory 

The first example is the quadratic f{R) gravity given in section (j2.4p . This is a particular 
simple example, since -F is a linear function with a simple inverse. Using the procedure 
above, we find that the gravity/scalar theory in Einstein frame is given by 

C = ^{R-\{dckf-V), (107) 
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where the scalar potential can be expressed in terms of a superpotential, namely 
We find that 

= "^ = ~ 2V2(Z^-''l)(I^-2) - ^'"'^ 

The domain wall solution in the f{R) gravity can be obtained by treating R as the 
coordinate, rather than z. In other words, we have 

where 

R + AD{D-l)W{Rl 

I(D^) ■ (^^^^ 

The function A is now given by 

f 2WWr , , , 

The domain wall approaches the AdS boundary at R = Rq where X{Ro) = 0. 
5.2 Another example 

Another example is provided by (j27p . In D = 4, we have 

£4 = R^R + P. (113) 
The resulting scalar /gravity system (|107p has a complicate scalar potential, given by 



V = -^^"^ [^^ -3(3+ v¥4 + 3^j - 3$ + Y 2$2 + 3/3 + ^/^^+3P^j , (114) 
where $ = e^'i'l^ . 

It should be pointed out that for the majority of our /(i?) gravities that admit Killing 
spinor equations, it is unnatural to convert them to the Brans-Dicke theory. If one insists 
on doing so, the philosophy should be applied to Einstein gravity with a Gauss-Bonnet 
term where the R + aF? part should be converted to the Brans-Dicke theory as well. The 
consequence is that Einstein gravity with Gauss-Bonnet term should be viewed as the Brans- 
Dicke theory coupled with the Ricci and Riemann tensor square terms. This formalism is 
clearly less elegant than the original pure gravity formalism. 
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6 Linear spectrum in (A)dS 

As was discussed in section 3, f{R) gravity admits (A)dS metrics as its vacuum solutions. 
We have constructed a large number of "BPS" domain wall and cosmological solutions that 
run from one (A)dS to the other. It is a formidable task to examine the stability of these 
solutions. In this section, we study the linear fluctuation of f{R) gravity in such a (A)dS 
vacuum instead. As has been discussed in section 2, there are two types of (A)dS vacua that 
could arise in f{R) gravities. The first type is the usual one that satisfies ([5]). The second 
type is the one we discovered in this paper and it is characterized by the divergent F{Rq). 
In this section, we shall be only concerned with the linearization f{R) gravities around the 
(A)dS vacua of the first type. Linearized f{R) gravity in such AdS4 were studied in [42] . 
For the linear perturbation g^u g^iu + h^y, we impose the gauge condition 

V^h^y = Vuh. (115) 

This gauge condition is different from the usual de Bonder gauge, but it is more effective 
to use in theories with a cosmological constant since it implies the vanishing of the trace 
scalar mode in Einstein gravity with a cosmological constant. It has been adopted in recent 
studies in critical gravities [43 ^ 1441 [45] . We find that the linearized equation of ([2]) becomes 

2A 



-\f{R^)[U-^^)H^, = 0, (116) 
-{D-l)kf"{Ro)iu-m')h = 0, (117) 



where 



n = h -L . 2(D - l)f"{Ro) ^ 



2_ i?o/"(i?o)-i(^-2)/'(i?oJ 



To derive the above equations, we have made use of the following formulae 



[D,V^]h = AV,h, [□,v^V,]/i = -^^J^,. (119) 



It is clear that H^i, is traceless; it is also transverse by the virtue of the equation of motion 
for h. In the above, we assume that Rq, /(Rq), f (Rq) and /"{Rq) are all non-vanishing. 
Thus we see that in general, in addition to the massless spin-2 graviton mode, there is also 
a massless scalar trace mode. However, there is no higher-order propagator for both modes, 
unlike the case in theories with more general higher curvature invariants. This is consistent 
with the fact that in terms of physical degrees of freedom, f{R) gravity is equivalent to 
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a special class of the Brans-Dicke theory. The lacking of higher-order propagators implies 
that there is no critical phenomenon as those discussed in [431 H5] . 

In the special case, where /"{Rq) = 0, the equation pi7p implies that h = 0, and the 
theory contains only the massless graviton, as in the case of Einstein gravity. If f'{Ro) = 0, 
graviton H^^, no longer has its kinetic term. If f'{Ro) = = f"{Ro), the theory has no 
propagating mode at all. For example, the theory 

^3 



C4 = V^{R-Ror (120) 

satisfies the criteria. What is interesting is that although such a theory does not have any 
perturbative propagating degrees of freedom, it nevertheless admits the (A)dS Schwarzschild 
black hole solution. This particular aspect of the theory is similar to three-dimensional 
Einstein gravity with a cosmological constant. 

For general case with non-vanishing f'{Ro) and /"{Rq), the ghost-free conditions are 

/'(i?o)>0, Rof{Ro)>0. (121) 

The tachyon-free Breitenlohner-Freedman (BF) condition in AdS is given by 

> ^Ro , (122) 

We now examine the stability of some of our f{R) gravities in the (A)dS vacua. The 
first example to consider is the quadratic Ricci-scalar action (124p . We have demonstrated 
in section 5 that this is equivalent to (jl07p . We choose a convention that the AdS fixed 
point for the superpotential ()109p occurs at = 0, which implies that 

{D-4){5D-2) 
^ = 3D{D + 2) • ^^2^) 

Expanding the scalar potential V around (p = 0, we find that 

V = ^^^Ro + iMV' + • • • , (124) 

where 

" 3D{D + 2)' l6{D-iy ^ ' 

It follows from (j26p that the reality condition requires that a < 0, and hence the vacuum 
is AdS. Furthermore, we have 

^ - - -16D(I)-1)3^° > °' 
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hence the BF bound is satisfied. Thus we find that the scalar /gravity theory is both tachyon 
and ghost free in the AdS vacuum with R = Rq. Now let us examine the corresponding 
f{R) gravity. For D ^ 4, there are two AdS vacua, namely 

R = Ro and R = Ro= ^.^~^. , (127) 

6{L> + 2)a 

where Rq is given by (jl25p . The R = Rq vacuum is "BPS", whilst the R = Rq one is not. It 
is easy to verify that F{Rq) = 1 and Rof"{Ro) = IuRq > 0. The calculated from (jllSp 
is exactly the same as M^. The situation is quite different for the non-"BPS" vacuum with 
R = Rq. Although there is no tachyon, the spin-2 graviton is a ghost field. If we reverse the 
overall sign of the action, the spin-0 trace mode becomes a ghost. That the "BPS" vacuum 
is stable whilst the non-"BPS" vacuum is unstable is consistent with our expectation. 

The second example we would like to examine is the two theories given in section 2.2. 
Both theories ([6]) and (jlip can give rise to the same cosmological solution ([8]) which describes 
an evolution from the flat spacetime to the inflationary de Sitter vacuum. It is clear that in 
the flat region, ([6]) is a good perturbative theory. However, in the region where R = Rq = 
48/3^, the theory becomes singular with divergent F{Rq) and /"{Rq). On the other hand, 
theory (jlip is opposite. In the R = region, the theory is singular, but it is well behaved 
in the R = Rq region, with F{Rq) and /"{Rq) given in (112p . and hence the vacuum is ghost 
free provided that /3(T2 > 0. Note that the same conclusion also holds for the pair of theories 
(|73]1 and ([52]) in D = 4. Furthermore the BF bound for the theory ([73]) in the "BPS" AdS 
vacuum is satisfied, and hence the theory is both ghost and tachyon free. 

Thus we see an interesting phenomenon in our f{R) gravities, which we have mentioned 
in section 2. For a giving Killing spinor equation, and hence one "BPS" domain wall or 
cosmological solution, there can be two f{R) gravities. For a solution that connects two 
different AdS vacua with A+ and A_, one f{R) theory is well-defined in the A-|_ vacuum 
with no ghost and tachyon, but becomes singular at A_, and vice versa for the other f{R). 

The third example we consider is the five-dimensional f{R) theory ()70p that admits the 
AdS wormhole solution ()69p . The AdS wormhole is symmetric with both AdS boundaries 
having the same Rq = —500. It is clear that the theory associated with o"2 is singular. On 
other hand, the theory associated with ai, namely 

f{R) = aiV-{R + 200)^, (128) 

is well defined in the vacuum. We have F{Rq) = IbVSai and Rof"{Ro) = 25^/3cri/2, 
and hence the vacuum fluctuation is ghost free for positive ai. The mass square of the 
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spin-0 mode is given by m? = —100, which precisely saturates the BF bound. This theory 
is different from the second example, in that both AdS boundaries have the same Rq and 
hence one theory is needed instead of having to have both theories to patch different regions. 

The last example we shah examine is the smooth Randall-Sundrum II solutions discussed 
in section 3. For the simpler case (j62p in D = 5, both F{Rq) and Rq/'^Rq) are given in 
(j66p . They cannot be both positive and hence the AdS vacuum suffers from having a ghost 
field. (It is easy to obtain m? = 300 for the scalar mode, and hence it is not a tachyon.) 
Of course, this is only one example of many possible RS II solutions, and it is of interest to 
investigate whether such a ghost problem of the "BPS" RS II in f{R) gravities is generic 
or not. Furthermore, we have imposed that f"{Ro) be finite, which is not entirely clear to 
be necessary. 

7 Conclusion 

In this paper, we follow the procedure outlined in [2011211122] and obtain the condition on the 
subclass of f{R) theories that admit Killing spinor equations. We present many examples 
of such f{R) gravities. One advantage of our theories is that the Killing spinor equations 
reduce the fourth-order Einstein equations for the domain wall and FLRW ansatze to very 
simple first-order equations, and hence exact solutions can be constructed. 

For domain wall solutions, we find exact smooth examples that describe the RS II sce- 
nario, AdS wormholes and the RG flow from the IR to the UV. In all these solutions, the 
metric runs from one AdS to another. This is very different from other higher-derivative the- 
ories such as Lovelock gravities with the Gauss-Bonnet term, which also have multiple AdS 
vacua, but have no known flow running from one to the other. Our examples demonstrate 
that f{R) is a fruitful arena to investigate and apply the AdS/CFT correspondence. 

Rich classes of exact and smooth cosmological solutions also emerge in our f{R) gravi- 
ties. We flnd a semi-realistic cosmological solution that evolves from an inflationary starting 
point to end with an ever-lasting expanding universe with a much smaller cosmological con- 
stant. We also flnd a pre-big bang model where a flat universe bursts into inflation by the 
non-perturbative effect of the higher-order curvature terms. In addition, we flnd smooth 
bouncing and crunching universes. 

Since the cosmological evolution in our f{R) gravities is solely governed by the equation 
(|84p. it is a matter of flnding the right proflle of W{R) in order to flt the observational 
data. However, one technical drawback is that for a given W, the f{R) is not determined 
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directly, but via a second-order linear differential equation, which may not have a close- form 
solution. Nevertheless we have obtained many explicit examples in this paper. Classically, 
it can be argued that this is not essential since W gives all the information. At the quantum 
level, the exact form of f{R) is likely to become much more important. It is of interest to 
investigate whether it is possible to compute the quantum effect on the information given 
by the W alone. 

The full analysis of the stability of our "BPS" solutions is beyond the scope of this 
paper. Instead, we investigate the stability of the (A)dS vacua that these metrics connect 
to. We study the linearized gravity around the (A)dS vacua. We adopt the gauge that was 
used previously for studying critical gravities [431 144j . In general, f{R) gravity consists of 
one massless spin-2 and one-massive spin-0 modes. We also obtain the condition for which 
the spin-0 mode decouples so that the spectrum is identical to that of Einstein gravity. 
More exotic situation can arise where a theory has no propagating degree of freedom, yet it 
admits the Schwarzschild (A)dS black hole as a solution, analogous to Einstein gravity with 
a cosmological constant in three dimensions. We obtain the conditions for f{R) theories to 
be absent from the ghost and tachyon fields, and give a detail analysis for a few examples. 

There is an intriguing phenomenon in our f{R) gravities. As we have mentioned, for a 
giving W in the Killing spinor equations, and hence one "BPS" domain wall or cosmological 
solution, there can be two f{R) gravities. For the solution that connects two different AdS 
vacua with A+ and A_, We find examples that one f{R) theory is well-defined in the A_|_ 
vacuum with no ghost and tachyon, but becomes singular at A_, and vice versa for the other 
f{R)- This suggests that there can exist multiple classical f{R) gravities that give the same 
full cosmological evolution; however, different stages of the evolution may select different 
specific theories for the quantum description. This is similar to the common phenomenon 
in differential geometry that a typical manifold requires multiple different but overlapping 
coordinate patches in order to cover it. 

It was shown in [461 147j that non-supersymmetric theories that admit Killing spinor equa- 
tions can be pseudo-supersymmetrized by introducing pseudo fermionic partners. In these 
theories, it can be shown that the Lagrangian is invariant under the pseudo-supersymmetric 
transformation rules up to the quadratic order in fermions. This suggests that there should 
be pseudo-supersymmetric versions of our f{R) gravities. It is of great interest to construct 
such f{R) pseudo-supergravities. 

To conclude, our construction of f{R) gravities that admit Killing spinor equations 
allows us to find exact "BPS" domain wall and FLRW cosmological solutions with varying 
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Ricci scalar R. The significance of these solutions is that they explore the function f{R) in 
contrast to the previously known solutions with fixed R. This opens a new door to study 
both the AdS/CFT correspondence and cosmology in the context of f{R) gravities. Our 
construction is based on f{R) theories in the metric formalism, and hence it is natural to 
extend our discussion to the Palatini formalism where both the metric and the connection 
are assumed to be independent variables. It is also of interest to investigate whether the 
Killing spinor equations of our f{R) gravities can be extended to include matter. 
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